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The  flow  result!):^  from  the  passage  of  a  shock  wave  through  a  dusty- 
adr  layer  Is  stvidled  analytically  and  numerically.  For  the  cases  treated 
here,  a  shock  wave  Is  always  reflected  at  the  first  contact  surface 
separating  the  pure  gas  from  the  dtisty  layer  and  a  shock  wave  is  trans¬ 
mitted  into  the  dusty  layer.  The  transmitted  shock  wave  is  stronger. 

The  criteria  for  the  reflected  wave  in  terms  of  the  properties  of  the 
air  and  dusty  layer  are  obtained  based  on  an  idealized  equilibrium^gas 
approximation.  Working  curves  are  presented  showing  the  effects  of  the 
svispended  particles  on  the  resulting  flow.  Similarly,  at  the  second 
contact  front  of  the  dvisty  layer,  a  rarefeuition  wave  is  reflected  and 
the  transmitted  shock  wave  transMtted  into  the  air  is  weakened  by  this 
nonlinear  interaction.  The  rarefaction  wave  reflects  at  the  first  con¬ 
tact  front  as  a  cospression  wave  and  proceeds  through  the  layer  to 
refract  at  the  second  contact  surface  where  it  reflects  as  a  rarefaction 
wave  and  a  transmitted  conpression  wave.-_Thls  wave  overtakes  the  trans¬ 
mitted  shock  wave  in  air  and  produces  th^final  transmitted  wave,  a  new 
contact  surface  and  a  weak  reflected  rar«fbction  wave.  This  final 
emergent  shock  wave  from  the  dusty  air  has  almost  the  same  strength  as 
the  original  shock  wave  entering  the  layer.  A  particular  case  for  an 
initial  shock  wave  pressure  ratio  of  5  has  been  chosen  to  illustrate 
this  type  of  interaction  in  detail. 

Finally,  the  time-dependent  transition  properties  through  the  shock 
waves,  contact  svirfaces  and  rarefaction  waves  were  found  by  solving  the 
equations  of  motion  numerically  using  a  modified  random-choice  method 
with  an  operator-splitting  technique.  This  provides  the  details  for  the 
formation  of  the  idealized  equilibrium  flows  in  the  dusty-air  layer  and 
the  reflected  and  transmitted  shock  waves  in  the  air  surrounding  the 
dusty-air  layer. 
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equilibrium  speed  of  sound 
C  specific  heat  of  solid  material 

SI 

gas  specific  heat  at  constwt  pressure 
Cy  gas  specific  heat  at  constant  volume 

specific  heat  at  constant  volume  of  idealized  equilibrium  gas 
D  drag  force  acting  on  a  particle 

d  particle  diameter 

k  gas  thermal  condiactivity 

m  mass  of  a  particle 

p  gas  pressure 

Q  rate  of  heat  transfer  to  a  particle 

R  gas  constant 

T  gas  teuqjerature 

t  time  from  shock  penetration  across  a  contact  surface 

u  gas  velocity 

u  propagation  velocity  of  shock  wave 
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V  velocity  of  particles 

X  space  coordinate  measured  from  a  contact  siirface 

6  shock  transition  thickness 
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drag  coefficient  of  a  particle 
shock  Mach  number 
Nusselt  nuniber  of  a  particle 
normalized  pressure 
gas  Prandtl  number 
particle  Reynolds  number 
normalized  teiq>erature 
normalized  velocity 
normalized  space  coordinate 
mass  concentration  ratio  cr/p 
specific  heat  ratio  of  two  phases 
gas  specific  heat  ratio 

specific  heat  ratio  of  idealized  equilibrium  gas 
normalized  mass  concentration 
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1.  mrooDUCTioM 


High-speed  flews  of  a  nlxtirre  of  a  gas  and  snail  solid  particles 
are  encountered  in  several  branches  of  engineering  and  science  (Refs. 

1-6) .  Shock  waves  in  such  a  dusty  gas  exhibit  relaxation  feattires 
typical  of  two-phase  flow.  For  sufficiently  strong  shock  waves,  a 
discontinuous  Juitp  in  the  gas  phase  precedes  a  thick  transition  region. 

A  nonequilibri\an  state  of  the  mixture  gives  way  gradually  to  an  equili¬ 
brium  tramsition  through  the  transfer  of  momentum  suid  heat  between  the 
gas  and  the  particles.  This  feature  was  utilized  in  ejfperimental  studies, 
using  a  shock  tube,  for  an  understanding  of  the  basic  mechanisms  of  the 
interaction  between  the  two  phases  (Refs.  7-9) •  For  example,  effective 
drag  coefficients  of  the  particles  were  deduced  from  the  experimental 
results, but  their  values  cannot  as  yet  be  considered  as  definitive. 

Some  experiments  used  a  shock  tube  in  which  only  a  limited  portion 
of  the  channel  was  filled  with  a  dusty  gas  (Refs.  10-l4).  Apparently 
a  uniform  dusty  gas  is  easier  to  produce  in  such  a  configuration.  It 
also  has  the  advantage  that  disturbances  arising  from  diaphragm  rtpture 
are  avoided.  In  this  case,  a  shock  wave  initially  propagating  in  a 
pure  gas  then  penetrates  the  dusty-gas  layer  and  gives  rise  to  a  refrac¬ 
tion  problem, 

A  nimierical  analysis  of  this  type  of  flow  was  done  by  Outa,  Tajima 
and  Morii  (Ref,  12),  using  the  method  of  characteristics.  However,  they 
made  an  error  in  calculating  the  dusty-gas  flow  by  assuming  that  the  mass 
concentration  of  the  particles  vanishes  along  the  boundary  psirticle  path. 
There  must  be  a  discontinixity  in  the  mass  concentration  at  the  boundary 
of  the  particles,  since  diffusion  was  not  considered.  Recently,  Marconi, 
Rudman  and  Calia  (Ref.  15)  solved  a  problem  of  shock  induced  flow  in  a 
dusty-gas  layer  using  a  characteristic-based  finite-difference  method. 

This  was  a  less  practical  case  since  the  incident  shock  wave  was  produced 
by  a  moving  piston  and  a  subsequent  reflected  shock  wave  at  the  piston 
affected  the  flow. 

In  the  present  analysis  the  basic  features  are  studied  of  the  flow 
induced  by  a  shock  wave  as  it  penetrates  into  a  dusty-air  layer.  The 
transient -flow  behaviour  is  clarified  by  solving  the  equations  of  motion 
using  a  modified  random-choice  method  (Ref.  l6)  together  with  an  operator¬ 
splitting  technique.  Discussions  based  upon  the  idealized- equilibrium-gas 
approximation  are  also  made  for  a  more  conplete  understanding  of  the  flow 
characteristics.  The  inverse  problem  is  also  considered  as  the  shock  wave 
leaves  the  dusty-air  layer  again  into  free  air.  The  strengths  of  the 
initial  shock  wave  as  it  penetrates  the  layer  and  the  emergent  shock  wave 
as  it  leaves  the  layer  are  conpared  and  discussed.  It  is  shown  that  the 
analysis  reduces  to  the  case  of  a  gaseous  layer  embedded  in  air  in  a 
shock- tube  channel,  very  similar  to  the  case  studied  by  Glass  and  co¬ 
workers  (Refs.  17-20)  nearly  three  decades  ago.  The  case  of  a  shock  wave 
passing  throtigh  a  helium  layer  embedded  in  a  channel  containing  air  is 
illustrated  in  the  time-distance  plane  on  Plate  1  and  its  explanatory 
sketch.  In  this  case  a  rarefaction  wave  is  reflected  at  the  first  contact 
sixrface  and  a  shock  wave  at  the  second  contact  surface,  xxnllke  the  passage 
through  a  dusty  layer  where  the  reverse  happens.  Hewi  argon  or  carbon-dioxide 
been  used,  then  it  would  hawe  been  similar  to  the  dusty-gas  layer. 
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Finally,  it  is  worth  noting  that  in  the  present  analysis  the  over¬ 
taking  of  the  emergent  shock  wave  by  the  refracted  cotqpresslon  wave  is 
also  considered  in  order  to  give  a  definite  answer  as  to  whether  or  not 
the  emergent  shock  wave  is  weaker  than  the  initial  shock  wave  entering 
the  dusty-air  layer.  In  the  particular  illustrative  exaag>les  treated 
here  for  an  incident  shock  pressure  ratio  of  5,  the  final  transmitted 
shock  wave  after  the  nonlinear  wave  interactions  ultimately  comes  out 
with  almost  the  same  strength  as  the  incident  one.  Of  course,  other 
exaiig>les  may  be  treated  by  the  reader  following  the  methods  outlined  in 
this  report. 


2.  I»RACTION  OF  A  SHOCK  WAVE  AI  AN  AIR/DUSTY-AIR  CONTACT  FRONT  FOR 
EftUniatlUM  FLOWS 

In  this  section,  we  consider  the  variation  of  a  shock  wave  when  it 
passes  acTOBS  the  contact  surface  separating  a  pure  gas  and  a  dusty  gas 
of  seml-lnfinlte  length.  This  corre^onds  to  the  case  when  the  thickness 
of  the  dusty-gas  layer  is  much  longer  than  a  characteristic  relaxation 
length.  As  time  elapses  after  the  shock  wave  passes  across  the  contact 
sxirface,  the  particles  become  in  equilibrium  with  the  gas  except  for  the 
transition  region  of  the  shock  wave.  The  dusty  gas  can  be  treated  effec¬ 
tively  as  a  perfect  gas  if  we  neglect,  from  a  large-scale  point  of  view, 
this  nonequilibrium  region  together  with  a  contact  region  of  finite  thick¬ 
ness,  which  may  arise  as  in  the  case  of  a  dusty-gas  shock  tube  (Ref,  l6) . 
The  specific  heat  at  consteint  volume,  the  specific  heat  ratio  and  the 
speed  of  sound  of  the  idealized  equilibrium  gas  are  given  by  (Ref.  4) 
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A  shock  wave  or  a  rarefaction  wave  is  reflected  at  the  contewt  sxurface 
separating  perfect  gases  of  different  kinds  when  an  incident  shock  wave  is 
transmitted  (Ref.  19,  see  Fig.  1),  The  type  of  reflected  wave,  i.e., 
shock  wave,  Mach  wave  or  rarefaction  wave,  is  determined  according  to  the 
conditions  (Ref.  17), 

respectively.  This  is  identical  to  P2/P2|  ^  i* 

First,  we  consider  a  shock  wave  transmitted  from  a  pure  gas  to  a  dusty 
gas.  For  this  case,  we  have  (see  Fig.  l), 
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If  we  assume  that  the  tenperatiares  are  the  scune  in  both  gases,  the  left- 
hand  side  of  Eq.  (^t)  is  calciilated  as  follows: 
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(6) 


Three  cases  can  be  found  from  the  above  expression  according  to  the 

value  of  p: 

(i)  For  P  ^  7/(7-!) >  the  right-hand  side  of  Eq,  (6)  is  positive  inde¬ 
pendent  of  the  values  of  a  and  P4/P5.  Therefore,  the  reflected 
wave  is  a  shock  wave. 

(ii)  For  (7+l)/(7-l)  >  P  >  7/(7-!),  two  cases  exist  depending  on  the 
value  of  g: 

(a)  For  g  ^  7-I-7/P,  the  last  factor  of  Eq.  (6)  is  positive 
because  P4/P5  is  greater  than  unity.  Therefore,  the  reflected 
wave  is  a  shock  wave  independent  of  the  strength  of  the  inci¬ 
dent  shock  wave. 

(b)  For  g  <  7-1-7/p,  there  exists  a  critical  pressure  ratio 
PU/P5  =  (P-l)[(2cep/7-l)  +  (7+1/7-1)  -  P]"l  for  which  the 
last  factor  of  Eq.  (6)  vanishes.  If  the  incident  shock  wave 
is  stronger  than  this  critical  one,  the  reflected  wave  is  a 
shock  wave.  On  the  other  hand,  a  rarefaction  wave  is  re¬ 
flected  if  the  incident  shock  wave  is  weaker  than  the 
critical  one. 

(iii)  For  P  >  (7+l)/(7-l),  three  cases  exist  d^ending  on  the  value  of 

a: 


(a)  For  g  >  7-1-7/p,  the  reflected  wave  is  a  shock  wave  independent 
of  the  strength  of  the  incident  shock  wave  as  in  (iia) . 

(b)  For  7-1-7/p  >a  >  (7-l)/2-(7+l)/2p,  a  critical  pressure  exists 
as  in  (iib).  A  shock  wave  or  a  rarefaction  wave  is  reflected 
according  to  the  strength  of  the  ixxjident  shock  wave. 

(c)  For  g  <  (7-1)/2-(7+1)/^,  the  last  factor  of  Eq.  (6)  is  negative 
independent  of  P4/P5.  Therefore,  the  reflected  wave  is  a  reire- 
f action  wave  for  any  incident  shock  wave. 
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When  a  shock  wave  is  transmitted  into  the  dusty  gas,  the  gas  is 
decelerated  by  the  particles  of  large  inertia.  This  deceleration 
tends  to  raise  the  pressvire.  On  the  other  hand,  the  particles  absorb 
heat  from  the  gas  at  high  teii5)erature.  This  absorption  of  heat  or 
cooling  of  the  gas  reduces  the  pressvure.  Thus,  if  the  specific  heat 
of  the  particles  is  so  small  that  the  absorbed  heat  is  sufficiently 
small,  the  pressxire  of  the  dusty  gas  becomes  higher  than  that  of  the 
pure  gas.  As  a  result,  a  reflected  shock  wave  is  induced,  for  example, 
as  for  case  (i).  Conversely,  for  particles  of  laurge  specific  heat,  a 
rarefaction  wave  appears  as  for  case  (iiic) . 

Generally,  the  specific  heat  ratio  of  the  two  phases  P  is  not  so 
high  and  case  (i)  becomes  practical,  that  is,  a  shock  wave  is  reflected 
at  the  contact  surface  in  most  cases.  However,  case  (iii)  can  exist. 
For  exanple,  consider  a  mixtiare  of  amgon  and  polystyrene-particles, 
for  which  y  =  1.6?  and  p  =4.4.  For  this  mixture,  we  can  have  the 
three  cases;  (iiia)  for  a  ^  0.29,  (iiib)  0.29  >  o:  ^  0.032,  and 
(iiic)  for  a  <  0.032.  Although  the  case  of  a  reflected  rarefaction 
wave  is  interesting,  we  have  found  that  the  rarefaction  wave  is  very 
weeik  for  this  mixture  and  it  will  not  be  discussed  further . 


For  the  case  of  a  reflected  shock  wave,  the  pressure  behind  the 
wave  P3(=P2)  is  given  by  (Refs.  17,  19) 


Several  working  ciorves  for  the  mixt\xre  of  air  (7  =  1.4)  and  glass 
particles  (P  =  l)  are  presented  illustrating  the  effect  of  the  exis¬ 
tence  of  particles  upon  the  uniform  states  between  the  transmitted 
shock  wave,  the  contact  surface  and  the  reflected  shock  wave.  Figure 
2  shows  the  variations  of  P3/P5(=P2/Pl)  with  the  pressure  ratio  of 
the  incident  shock  wave  P4/P5  for  values  of  a  over  the  range  0  <  a  <  2. 
The  pressure  po  increases  with  a  for  a  fixed  P4/P5,  that  is,  the  “ 
reflected  shock  wave  and  transmitted  shock  wave  are  stronger  for  larger 
mass  concentration  of  particles. 

Other  flow  quantities  are  calculated  from  the  Rankine-Hugoniot 
relations  (see  Appendix  A).  The  variations  of  the  velocity  of  the 
fluid  between  the  reflected  and  transmitted  shock  waves  "^2/^  with 
the  incident  shock  pressure  ratio  are  shown  in  Fig.  3.  The  particles 
decelerate  the  gas  and  therefore  the  mlxtiire  velocity  reduces  with  a 
for  a  fixed  P4/P5.  The  ten^erature  behind  the  transmitted  shock  wave 
T2  decreases  as  a  result  of  absorption  of  energy  by  the  particles 
(Fig.  4),  but  the  pure  gas  is  heated  by  the _ reflected  shock  wave 


(Fig.  5).  The  density  ratios  P2/P5(^a^o’l)  P3/P5  ®*’®  raised  by 

the  deceleration  of  the  gas  due  to  the  existence  of  the  particles 
and  by  the  induced  reflected  shock  wave,  respectively  (Figs.  6,  7). 

The  former  is  larger  than  the  latter,  as  expected.  A  contact  region 
s^arates  these  layers  of  different  d^sities  and  ten^eratures. 

Figure  8  shows  the  variations  of  the  transmitted  shock  Mach 
number  based  on  the  frozen  speed  of  sound  equal  to  a^.  The  presence 
of  the  particles  reduces  the  speed  of  propagation  of  the  transmitted 
shock  wave.  For  sufficiently  high  mass  concentration  of  the  particles, 
the  transmitted  shock  wave  propagates  at  a  lower  velocity  than  the 
frozen  speed  of  sound  so  that  it  must  be  fully  dispersed. 


3.  REFRACTION  OF  A  SHOCK  WAVE  AT  A  PIETY- AIR/AIR  COMTACT  PROMT  FOR 
EQ.UILIBRIUM  FLOWS 

Next,  we  consider  a  shock  wave  transmitted  from  a  dusty  gas  to  a 
pure  gas.  For  this  case,  we  have 


II 

Cvl  =  "• 

" 

(8) 


The  left-hand  side  of  Eq.  (4)  is  calculated  \ri.th  these  values  in  a 
similar  way  to  Eq.  (6) .  The  same  factor  as  the  last  one  on  the  right- 
hand  side  of  Eq.  (6)  also  appears,  but  with  a  changed  sign.  Therefore, 
the  criteria  (i),  (ii)  and  (iii)  explained  above  still  holds  for  this 
case,  although  the  reflected  shock  wave  should  now  read  as  the  reflected 
rarefaction  wave  and  vice  versa.  Thus,  a  rarefaction  wave  is  reflected 
generally  at  the  contact  surface. 

When  a  rarefaction  wave  is  reflected  at  the  contact  surface,  the 
pressure  behind  the  wave  P3(=P2)  is  given  by  (Refs.  17,  19) 


We  now  present  the  working  curves  showing  the  variations  of  the 
flow  quantities  with  the  incident  shock  pressure  ratio  P4/P5  for  the 
mixture  of  air  and  glass-particles.  The  pressure  behind  the  trans¬ 
mitted  shock  wave  P2  is  smaller  than  pi4  behind  the  incident  shock 
wave  in  the  dusty  gas  (see  Fig.  9)  and  a  rarefaction  wave  is  reflected. 
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Variations  of  the  velocity  In  the  region  between  the  transnltted 
and  the  reflected  waves  nomslized  by  the  frozen  speed  of  soiind 
equal  to  a^  are  shown  in  Fig,  10.  The  values  of  the  velocity  ui,  of 
the  dusty  gas  behind  the  incident  shock  wave  are  also  shown  by  the 
thin  lines  for  coo5>arison.  The  gas  in  state  (3)  is  accelerated  to 
a  higher  velocity  stfter  it  is  released  from  the  restraints  of  the 
dust  particles.  The  tenperature  of  the  gas  behind  the  transmitted 
shock  wave  becomes  higher  than  that  behind  the  incident  shock  wave 
(Fig,  ll).  On  the  other  hand,  the  dusty  gas  cools  after  the  reflected 
rarefaction  wave  passes  through  (Fig.  12).  Larger  changes  in  ten?)er- 
ature  Eore  produced  for  larger  mass  concentrations  of  particles.  The 
densities  of  the  gases  behind  the  transmitted  and  the  reflected  waves 
are  smaller  than  that  behind  the  incident  shock  wave  (Figs.  13,  lU) . 
Comparing  these  two  figures,  we  see  that  the  gas  behind  the  reflected 
rarefaction  wave  is  denser  than  the  gas  behind  the  transmitted  shock 
wave. 


Variations  of  the  transmitted  shock  Mach  number  with  the  incident 
shock  pressure  ratio  are  shown  in  Fig.  15.  It  is  seen  that  the  trans¬ 
mitted  shock  wave  propagates  at  higher  speed  than  the  propagation 
velocity  of  the  incident  shock  wave  in  the  dusty  gas.  It  is  seen 
that  there  is  a  region  where  a  fully  dispersed  shock  wave  in  the 
dusty  gas  can  eventiaally  became  a  sharp  fronted  transition  in  air. 

That  is  the  reverse  of  the  case  of  refraction  from  air  into  a  dusty 
gas. 


4 .  TIME-DEPENDEMr  UOMEQUILIBRIUM  REFRACTION  AT  AN  AIR/DUSTY-AIR 
COMTACT  FRONT 


The  transient  flow  induced  after  a  shock  wave  passes  across  a 
contact  surface  is  discussed  in  this  section.  The  dusty  gas  is  assumed 
to  consist  of  a  perfect  gas  and  a  lot  of  solid  spheres  of  imiform  size. 

The  viscosity  and  thermal  conductivity  of  the  gas  are  taken  into  account 
only  for  the  interaction  between  the  gas  and  the  particles.  We  consider 
the  cases  when  the  mass  concentration  of  the  particles  is  conparable  to 
the  density  of  the  gas.  The  voliime  occupied  by  the  particles  is  negligible 
because  the  density  of  the  solid  particle  is  much  larger  than  that  of  the 
gas. 


The  equations  of  continuity,  momentum  and  energy  for  either  the  gas 
or  the  particles  are  given  by  (Refs.  1-6) 


?t 


If  (pu)  =  0 

(10) 

^  (av)  =  0 

(11) 

(12) 
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(13) 


If  (tv)  +  Ij  (tV^)  .2d 

|E{'’("n®*i''^)}+|:{Tv(c„etiv2)}=  2(^^)  (15, 

where  p,  p,  T,  u  are  the  pressure,  density,  tenperatiace  and  velocity 
of  the  gas  and  a,  0,  v  are  the  mass  concentration,  temperature  and 
velocity  of  the  particles,  respectively.  The  equation  of  state  for  the 
thermally  perfect  gas  is  given  by 


p  =  pRT 


(16) 


The  gas  and  the  particles  Interact  with  esush  other  through  the 
drag  force  D  and  the  heat  transfer  rate  Q  experienced  by  a  particle 
in  nonequilibrium  with  the  gas .  We  take  in  the  present  study  (Refs . 

21,  22) 

D  =  5  d^p(u-v)|u-v|Cjj 

=  J  d^p(u-v)  lu-v|(0.U8  +  28Re‘°'®^)  (17) 

Q  =  TrduC  Pr‘’^(T-6)Nu 

=  TrdpCpPr"\T-e)(2.0  +  0.6Pr%e^)  (l8) 

where  Re  is  the  Reynolds  number  based  upon  the  diameter  of  the  particle 
and  the  relative  velocity  of  the  particle  to  the  gas. 

Re  =  p(u-v|d4i  (19) 

and  Pr  is  the  Prandtl  number  of  the  gas, 

Pr  =»  ^Cp/k  (20) 

The  viscosity  and  the  thermal  conductivity  of  the  gas  vary  with 
the  tenperature.  We  consider  air  in  this  section,  for  which  (Ref.  23) 

p  *  1.71x10  X  f  ^  J  poise  (21) 

and 


Pr  »  0.75 


(22) 
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VIhile  the  equations  governing  the  motion  of  the  dusty  gas  have  been 
described,  the  pure  gas  obeys  Eqs.  (lO),  (l2)  and  (lU)  with  a  =  0, 
and  Eq.  (l6) . 

The  numerical  method  applied  to  the  present  problem  is  the  same 
as  in  Ref.  l6,  that  is,  a  modified  randonwchoice  method  together  with 
an  operator-splitting  technique  as  follows.  Two  sets  of  equations 
derived  from  Eqs.  (l0)-(l5)  are  considered.  One  set  is  obtained  by 
omitting  the  inhomogeneous  terms  on  the  right-hand  sides.  The  other 
set  is  obtained  by  omitting  the  terms  differentiated  with  respect  to 
X.  Solving  these  two  sets  of  equations  alternately  for  each  time  step 
results  in  the  solution  of  the  problem. 

The  former  set  of  homogeneous  equations  is  solved  by  the  rauidom- 
choice  method  (Refs.  l6,  24-27).  In  this  method  discontinuities  can 
be  described  definitively,  as  artificial  viscosity  is  not  required. 

The  elementary  calculation  of  the  random-choice  method  is  to  determine 
by  random  sampling  the  solutions  at  the  mid-point  between  two  adjacent 
points,  where  the  solution  is  known  at  a  previous  time.  The  equations 
for  the  gas  are  not  coupled  with  those  for  the  particles  in  this  set 
of  equations.  Taking  step-like  initial  conditions  for  the  gas  phase 
in  the  elementary  calculation,  one  obtains  classical  shock-tube  solu¬ 
tions  (Refs.  25-27).  On  the  other  hand,  one  assumes  initially  a 
linear  distribution  of  particles  in  order  to  avoid  the  difficulty 
that  the  solution  becomes  multivalued  for  step- like  conditions.  A 
st^-like  initial  condition  is  taken  only  at  a  boundary  of  the  particles 
for  the  boimdary  to  be  definite.  Further  details  are  given  in  Ref.  16. 

Numerical  calculations  were  done  for  the  mixture  of  air  (7  =  1.4) 
and  glass-particles  (P  =  l)  of  10  iim  diameter  with  a  mass  concentration 
ratio  a  =  1,  A  space  mesh  size  of  AX  ==  0.1  was  used  in  most  cases  and 
a  time  step  was  taken  as  the  Courant-Friedrichs-Lewy  condition  allowed. 

The  transition  of  the  flow  occurring  when  the  incident  shock  wave 
of  P4/P5  =  5  in  a  pinre  gas  penetrates  into  a  semi -infinite  dusty-gas 
region  is  shown  in  Figs.  16-19.  The  dusty  gas  was  assiuned  initially 
to  be  at  atmospheric  pressure  and  room  tenperature .  In  all  figures, 
the  distance  from  the  initial  contact  surface  is  measured  in  units  of 
4ppd/3pi  =  2.72  cm.  Thermodynamic  quantities  are  normalized  by  the 
corresponding  values  of  the  dusty  gas  before  the  shock  transmission. 

The  velocities  are  measured  in  the  speed  of  sound  for  the  pure  gas 
(7Pl/Pi)2  =  350  m/sec. 

Flow  structures  arising  when  a  small  time  (t  =  0.78x10“^  sec)  has 
elapsed  after  the  shock  penetrates  the  dusty  layer  are  shown  in  Fig.  16. 
The  particles  cannot  respond  quickly  enough  to  the  abrupt  change  in 
the  state  of  the  gas  at  the  shock  front,  but  are  gradually  accelerated 
and  heated  by  the  gas  around  them  (Figs.  l6c,  d).  On  the  other  hand,  the 
tenperature  and  velocity  of  the  gas  in  the  dusty-gas  region  are  reduced 
from  the  values  for  the  incident  shock  wave.  This  interaction  of  the 
two  phases  reduces  the  frozen  shock  juip  at  the  shock  front.  The  gas 
is  compressed  as  a  result  of  deceleration.  Figures  l6a  and  b  show  that 
the  pressure  and  density  of  the  gas  in  the  dusty-gas  region  are  raised 
higher  than  the  incident  values.  They  take  on  maxima  at  the  boiuidsury 


of  the  particles  at  this  stage.  The  conpression  of  the  gas  results 
in  a  shock  wave  that  is  reflected  back  into  the  pure  gas.  Its  initial 
fonoation  can  be  seen  in  Fig.  l6.  It  is  noted  in  Fig.  l6b  that  the 
particle  coiKsentration  rises  gradixally  behind  the  shock  front,  reeiches 
a  maximum  and  then  drops  suddenly  to  zero. 


Subsequent  transitional  behaviours  of  the  flow  are  shown  in  Fig. 

17  (t  =  1.56x10“^  sec)  and  Fig.  l8  (t  =  3 •12x10”*^  sec).  The  differ¬ 
ences  in  teniperature  and  velocity  between  the  gas  and  the  particles 
behind  the  transmitted  shock  wave  diminish  with  time.  While  the  dis¬ 
continuous  junp  in  the  gas  pheuse  at  the  transmitted  shock  front  decays, 
the  formation  of  the  reflected  shock  wave  in  the  pure  gas  becomes  more 
distinguished.  It  can  be  seen  from  Figs.  17b  and  c  that  a  contact  region 
including  the  boundary  of  theparticles  begins  to  form. 


The  variations  of  flow  quantities  at  t  ®  2.8lxl0"3  sec  are  shown 
in  Fig.  19.  The  exact  solution  for  the  idealized  equilibrium  gas  is 
also  presented  for  conparison.  The  present  results  agree  well  with  the 
idealized  solution.  The  reflected  shock  wave  in  the  pure  gas  is  almost 
steady,  although  it  is  still  followed  by  a  weak  conpressive  part.  Its 
position  is  a  little  behind  that  for  the  idealized  solution.  This 
difference  is  brotjght  about  by  a  delay  in  response  of  the  particles  when 
the  incident  shock  wave  penetrated  into  a  dusty  gas.  The  reflected 
shock  wave  moves  at  a  small  speed  in  the  laboratory  frame,  since  the 
induced  velocity  of  the  gas  bdiind  the  incident  shock  wave  is  fairly 
large,  A  stationary  structure  of  the  contact  region  is  seen  in  Figs. 

19b  and  c,  where  the  particles  are  in  equilibrium  with  the  gas.  The 
mass  concentrations  and  tenperatures  of  the  two  phases  vary  smoothly  in 
the  region  and  change  discontinuously  at  the  boiuidary  of  the  peu-ticles. 
The  structure  of  the  contact  region  reflects  the  history  of  the  inter¬ 
action  which  the  particles  experienced  so  far.  An  almost  stationary 
structure  of  the  transmitted  diock  wave  in  the  dusty  gas  has  developed, 
which  consists  of  an  initial  frozen  ^ock  front  followed  by  the  relaxa¬ 
tion  region.  The  particles  are  in  nonequilibrium  with  the  gas  in  the 
transition  region.  The  particles  and  the  gas  interact  with  each  other 
to  establish  an  equilibrium  profile. 

The  shock  pressure  ratio  P2/P1  is  found  to  be  6.09,  from  Eq.  (7), 
The  stationary  structure  of  the  transmitted  shock  wave  can  be  obtained 
separately  by  solving  the  set  of  ordinary  differential  equations  for 
a  coordinate  system  moving  steadily  with  the  shock  wave.  Figure  20 
shows  the  results  for  the  transition  region  obtained  by  making  use  of 
the  Runge-Kutta  method.  Conparing  Figs.  19  and  20,  we  see  the  structure 
of  the  transmitted  shock  wave  in  the  nonstationary  solution  is  close  to 
that  of  the  stationary  shock  wave.  It  should  be  noted  that  the  tenper- 
ature  of  the  gas  exhibits  an  wershoot  in  the  transition  region. 


Next,  we  consider  a  case  when  the  incident  shock  pressure  ratio  is 
so  small  that  a  fully  dispersed  shock  wave  may  arise  in  the  dusty-gas 
l^er.  Developments  of  the  flow  «rfter  the  incident  shock  wave  with 
Pl*/P5  =  1*5  penetrates  into  a  dusty-gas  region  of  send -infinite  length 
are  shown  in  Figs.  21,  22  and  23.  The  flow  structure  at  t  »  3 •12x10“^ 
sec  is  given  in  Fig,  21.  The  frozen  junp  at  the  shock  front  is  seen 
to  decrease  but  the  formation  of  the  reflected  shock  wave  in  the  piare 
gas  is  very  gradual  in  conparison  with  the  case  of  a  strong  shock  wave 
(see  Fig.  18). 
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The  variations  of  flow  quantities  at  t  =  2.8lxlO"3  sec  are  shown 
in  Fig.  22.  The  solution  for  the  idealized  eqxiilibrium  gas  is  also 
illustrated  by  the  dashed  lines  for  con^jari son .  The  present  results 
agree  fairly  well  with  it,  but  the  flow  is  still  developing.  The  flow 
quantities  vary  gradually  in  the  reflected  wave  and  a  discontinuous 
junp  does  not  form  as  yet.  The  frozen  front  of  the  transmitted  shock 
wave  has  almost  decayed. 

The  flow  structure  at  a  later  time,  t  =  5.30xl0“3  sec,  is  shown  in 
Fig.  23.  It  can  be  seen  that  the  coii5)ressive  reflected  wave  has  further 
steepened.  An  almost  stationary  contact  region  has  formed,  which  exhibits 
discontinuities  in  density  and  tenperature  at  the  boundary  of  the  particles 
It  should  be  noted  that  the  transmitted  shock  wave  has  dispersed  much  more 
since  t  *  2.8lxl0“3  sec  (note  that  differ  ant  scales  of  distance  are  used  in 
Figs.  22  and  23).  The  stationary  structure  of  the  fully  dispersed  shock 
wave,  for  which  the  shock  pressure  ratio  is  1.584  corresponding  to  the 
present  nonstationary  case,  was  solved  by  means  of  the  Runge-Kutta  method 
and  is  presented  in  Fig.  24.  A  comparison  of  the  shock  structures  in  Figs. 
23  and  24  shows  that  the  flow  over  the  transition  region  in  Fig.  23  does  not 
attain  the  final  stationary  state.  Further  dispersion  through  the  inter¬ 
action  between  the  gas  and  the  particles  must  follow.  Thus,  much  time  is 
needed  for  the  flow  to  become  stationary  even  after  the  frozen  shock  front 
has  decayed  ( see  Appendix  B) . 


5 .  TIME-DEPENDEIW  NONEQUILIBRIUM  REFRACTION  AT  A  DUSTY-AIR/AIR  CONTACT 

FRONT 

The  inverse  problem  will  be  considered  for  the  case  when  a  station¬ 
ary  shock  wave  in  a  dusty  gas  is  transmitted  into  a  pure  gas  of  semi-infinite 
length.  We  take  the  stationary  shock  transition  diown  in  Fig.  20  as  the 
initial  condition  for  the  incident  shock  wave.  The  flow  structure  at 
t  =  0,78x10"^  sec  after  the  shock  penetration  is  shown  in  Fig.  25.  The  gas 
ahead  of  the  particles  moves  at  a  higher  velocity  than  it  moved  in  the 
relaxation  region  of  the  incident  shock  wave,  since  the  gas  is  not  restrained 
by  the  particles  (Fig.  25d) .  This  increases  the  discontinuous  junp  at  the 
shock  front  transmitted  into  the  pure  gas.  Conciurrently,  an  expansion  of 
the  gas  is  induced  in  the  downstream  region  behind  the  discontinuous  shock 
front.  The  pressure,  density  and  teii5)erature  of  the  gas  decrease  from  the 
values  for  the  incident  shock  wave.  The  early  formation  of  the  reflected 
rarefaction  wave  can  be  seen  in  Fig.  25.  The  particles  are  overtaken  by 
the  gas  of  higher  velocity  and  tenperature  to  be  accelerated  and  heated 
(Figs.  25c,  d). 

The  flow  behaviour  at  t  =  1,56x10"**  and  3.12x10“**  sec  are  shown  in 
Figs,  26  and  27,  respectively.  The  discontinuous  transmitted  shock  wave 
increases  in  strength.  The  rarefaction  wave  reflected  into  the  dusty  gas 
continues  to  develop.  The  tenperature  and  velocity  of  the  particles  are 
increased  further.  It  should  be  noted  in  Fig.  27c  that  the  tenperature  of 
the  pairticles  is  higher  than  the  gas  tenperature  after  the  rarefaction  wave 
has  passed  owing  to  their  thermal  inertias.  Note  also  that  a  contact  region 
begins  to  form  ahead  of  the  boundary  of  the  particles  (Fig.  27b) . 

F'urther  develcpment  of  the  flow  is  shown  in  Fig.  28  for  t  =  2.8lxlO"3 
sec.  The  exact  solution  for  the  Idealized  equilibrium  gas  is  also  shown 
by  the  dashed  lines  for  conparison.  The  transmitted  wave  can  be  seen  to  be 


quite  stationary.  The  difference  in  its  position  between  the  present  result 
Eind  the  idealized  solution  is  due  to  the  delay  time  required  for  its  forma¬ 
tion.  It  is  also  seen  from  Figs.  28b  and  c  that  a  stationary  contact  region 
has  formed  around  the  boundary  of  the  particles.  This  contact  region  has 
no  discontinuity  in  the  density  and  tengjerature  of  the  gas.  The  reflected 
rarefaction  wave  in  the  dusty  gas  is  still  developing.  It  should  be  noted 
that  the  rarefaction  wave  is  carried  by  the  dusty  gas  flow  to  the  ri^t 
against  its  direction  of  propagation  owing  to  the  high  dust  velocity.  While 
its  position  agrees  well  with  that  of  the  idealized  solution,  the  rarefac¬ 
tion  wave  is  considerably  extended  by  conparison  \rith  the  idealized  solution 
owing  to  the  differences  in  ten5>erature  and  velocity  between  the  two  phases. 
The  calculations  were  done  vq)  to  t  =  5.9xlO“3  sec.  Although  at  later  time 
the  flow  structmre  exhibited  almost  similar  characteristics,  the  degree  of 
extension  of  the  rarefaction  wave  diminishes  relatively.  Since  the  head 
and  tail  of  the  rarefaction  wave  must  ultimately  propagate  at  the  equilibrium 
speed  of  sovind,  the  extent  of  the  wave  will  approach  that  of  the  idealized 
solution  after  a  long  time. 

The  transmission  of  a  fully  dispersed  shock  wave  from  a  dusty  gas  to 
a  pure  gas  is  solved  taking  the  stationary  shock  structure  shown  in  Fig.  24 
as  the  initial  condition  for  the  calculation.  The  flow  behaviour  at  t  = 
l.SyxlO'^  sec  after  the  shock  penetration  into  the  pure  gas  is  shown  in 
Fig.  29.  Only  small  chemges  in  pressure,  density  and  velocity  can  be  seen. 
The  tenqperature  of  the  gas  ahead  of  the  boundary  of  the  particles  is  raised 
and  a  sharp  peak  in  ten^ierature  is  observed  (Fig.  29c),  Figure  30  shows 
the  flow  structure  at  t  =  3«74xlO“3  sec.  The  flow  has  changed  considerably, 
but  is  still  far  from  the  idealized  equilibrium-gas  flew  illustrated  by  the 
dashed  lines.  The  variation  of  the  pressiye  does  not  exhibit  any  definite 
structvires  of  the  transmitted  and  reflected  waves.  The  temperature  and 
velocity  of  the  gas  ahead  of  the  particles  are  raised  to  the  values  higher 
than  the  incident  values.  However,  the  variations  of  flow  quantities  there 
are  very  gradual.  Steep  changes  in  density  and  temperature  near  the  boundary 
of  the  particles  indicate  the  formation  of  a  contact  region. 

The  flow  structure  at  a  later  time,  t  =  6.86xlO"3  sec,  is  shown  in  Fig. 
31.  Although  the  variation  of  the  flow  is  still  gradual,  we  can  distinguish 
the  formations  of  the  transmitted  shock  wave  and  the  reflected  rsirefaction 
wave.  The  form  of  the  transmitted  shock  wave  in  the  pure  gas  becomes  very 
steep.  Nonlinear  effects  will  make  this  wave  discontinuous  in  time.  The 
temperature  and  velocity  of  the  particles  have  approached  those  of  the  gas 
and  a  structure  of  the  reflected  rarefaction  wave  in  the  dusty  gas  can  be 
seen  in  the  figures.  The  width  of  the  rarefaction  wave  is  much  larger  than 
that  of  the  idealized  solution  as  in  the  case  of  a  strong  incident  shock 
wave.  The  structure  of  the  contact  region  esm  also  be  distinguished.  The 
uniform  states  of  the  flow  separated  by  the  tremsmitted  shock  wave,  the 
contact  region  and  the  reflected  wave  are  almcst  attained.  Their  positions, 
however,  are  much  different  from  those  of  the  idealized  equilibrium  gas. 

This  means  that  much  time  is  needed  for  the  relaxation  of  the  two  phases 
to  be  accomplished  when  the  change  of  the  flow  is  small. 


6.  NOHEftUILIBRIUM  PASSAGE  OF  A  SHOCK  WAVE  THROUGH  A  DUSTY-GAS  LAYER 

The  transmission  of  a  shock  wave  by  a  dusty  gas  layer  of  finite  thickness 
is  studied  next.  Figure  32  shows  the  flow  structure  at  t  =  2.18x10”^  sec 


after  the  shock  penetration  into  the  dusty-gas  layer  of  27,2  cm  initial 
thickness.  (The  dusty  gas  has  filled  initially  the  region  of  0  <  X  <  10 
in  the  figures.)  Other  conditions  were  the  same  as  those  of  the  semi¬ 
infinite  case  considered  in  Figs.  16-19.  It  is  seen  in  Fig.  32  that  the 
shock  wave  has  passed  through  the  layer.  The  thickness  of  the  layer  is  so 
thin  that  the  transmitted  shock  wave,  when  it  was  in  the  dusty-gas  layer, 
could  not  attain  its  stationary  value.  The  reflected  shock  wave  arising 
when  the  incident  shock  wave  passed  across  the  contact  surface  can  be  seen 
to  propagate  in  the  pure  gas  to  the  left  in  the  figures.  After  the  trans¬ 
mitted  shock  wave  passed  through  the  dusty-gas  layer,  a  rarefaction  wave 
is  induced  to  propagate  in  the  layer.  This  raorefaction  wave  penetrates 
into  the  pure  gas  upstream.  The  wave  is  seen  to  lie  between  the  reflected 
shock  wave  and  the  boundary  of  the  psurticles  in  Fig.  32.  Although  the 
rarefaction  wave  does  not  overtake  the  reflected  shock  wave  as  yet  in  the 
figures,  the  latter  will  almost  vanish  after  being  overtaken  by  the  former. 

The  magnitude  of  the  discontinuities  in  the  flow  variables  at  the 
frozen  shock  front  diminishes  while  the  shock  wave  propagates  in  the  dusty- 
gas  layeras  in  Figs.  l6-l8.  After  it  comes  out  of  the  dusty-gas  layer, 
its  magnitude  recovers  with  time.  It  can  be  seen  from  Fig.  32  that  the 
values  of  flow  quantities  immediately  bdiind  the  transmitted  shock  front 
are  very  close  to  those  behind  the  incident  shock  wave.  A  weak  con?)ression 
wave,  CW,  is  seen  to  follow  the  shock  front  in  Fig.  32a,  arising  when  the 
reflected  rarefaction  wave  passes  across  the  upstream  contact  surface  of 
the  dusty-gas  layer,  C  (Fig.  3^).  This  con5)ressive  wave  is  transmitted  out 
of  the  layer  and  overtakes  the  preceding  shock  front  Btg  .  Their  subsequent 
overtaking  cause  the  shock  wave  to  achieve  almost  its  original  incident 
strength. 

The  dusty-gas  layer  has  been  conpressed  to  be  9.2  cm  thick  as  is  seen 
from  Fig.  32b,  The  pressure  and  the  velocity  are  almost  uniform  over  the 
layer.  Only  the  mass  concentration  and  the  ten5)erature  vary  across  the 
layer,  but  the  particles  are  almost  in  equilibrium  with  the  gas  (see  Figs. 
32c  and  d) .  It  can  be  seen  from  Figs.  32b  and  c  that  the  nonuniform  state 
of  the  gas  protrudes  downstream  of  the  layer.  These  nonuniform  conditions 
of  the  flow  around  the  dusty-gas  layer  are  the  remains  from  the  passage  of 
the  shock  wave  through  the  layer.  The  nonimiform  state  will  last  until  the 
viscous  and  thermal  diffusion  of  the  gas,  as  well  as  the  particle  diffusion, 
remove  it. 

Figure  33  shows  the  flow  structure  at  t  =  i4,99xlO"3  sec  after  the  shock 
penetration  into  a  dusty-gas  layer  of  the  136  cm  initial  thickness.  (The 
initial  position  of  the  dusty  gas  was  over  0  <  X  <  50  in  Fig.  33) .  A  re¬ 
flected  rarefaction  wave  condng  out  of  the  dusty-gas  layer  is  clearly  seen 
to  follow  the  reflected  shock  wave.  The  transmitted  shock  wave  is  not 
affected  yet  by  a  compression  produced  when  the  reflected  rarefaction  wave 
passed  across  the  upstream  boundary  of  the  particles.  The  thickness  of 
the  dusty-gas  layer  is  fairly  large  so  that  the  downstream  flow  structure  is 
close  to  that  of  the  semi-infinite  case  (see  Fig.  28). 

When  the  thickness  of  the  dusty-gas  layer  is  very  large  (say  4  meters), 
we  will  find  the  reflected  shock  wave  in  the  pure  gas  is  like  that  shown  in 
Fig.  19.  The  reflected  reurefaction  wave  in  the  dusty-gas  layer  and  the 
transmitted  shock  wave  in  the  pure  gas,  on  the  other  hand,  will  be  as  shown 
in  Fig.  28,  Thus  the  strength  of  the  transmitted  shock  wave  (shock  pressure 
ratio  =  4,83)  is  weedcer  than  that  of  the  incident  shock  wave  (shock  pressure 


ratio  =  5.0)  if  we  neglect  the  overtaking  coiqpression  wave  that  can  result 
from  the  refraction  of  the  induced  rarefaction  wave  at  (see  Fig,  3^). 

The  transmitted  shock  wave  will  be  strengthened  when  it  is  overtaken  by  a 
shock  wave  developing  from  the  conpression  wave,  as  in  Fig.  3^.  The  shock 
pressxire  ratio  of  this  final  transmitted  shock  wave  can  be  calculated  on 
the  equilibrium- flow  approximation  and  is  found  to  be  5.005  for  the  present 
case.  Thus,  its  strength  for  all  intents  is  almost  the  same  as  that  of  the 
incident  shock  wave. 


7.  CONCLUSIONS 

The  passage  of  a  shock  wave  into  a  dusty-air  layer  was  solved  on  the 
basis  of  an  idealized  equilibrium  gas  which  ignores  the  transition  zones. 
Using  the  full  time-dependent  equations  of  motion,  it  was  shown  how  the 
shock  wave  and  contact  region  transitions  develop  to  their  equilibrium 
transitions.  The  idealized  equilibrium-gas  analysis  provided  the  criteria 
for  the  type  of  reflected  wave  which  occurs  at  the  contact  surface  separating 
a  pure  gas  from  a  dusty  gas.  In  the  present  case,  a  shock  wave  is  reflected 
when  the  incident  shock  wave  advances  from  air  into  a  dusty-air  layer  and 
a  rarefaction  wave  is  reflected  when  the  shock  wave  again  comes  out  of  the 
dusty  layer  into  air.  Working  curves  showing  the  effects  of  the  presence 
of  the  dusty-gas  layer  on  the  uniform  flow  states  formed  after  a  shock  wave 
passes  through  the  first  and  second  contact  surfaces  were  obtained  for  a 
dusty  gas  with  7  =  1.4  and  P  =  1.  The  strength  of  the  shock  wave  coming 
out  of  a  dusty-air  layer  of  large  thickness  can  be  found  from  these  curves 
as  a  function  of  the  incident  shock-pressure  ratio.  The  effects  of  the 
subsequent  wave  interaction  by  the  overtaking  conpression  wave  was  also 
analysed  for  one  particular  case  for  illustrative  purposes. 

The  transitions  of  the  flow  are  studied  numerically  by  means  of  a 
modified  random-choice  method  with  an  operator-splitting  technique.  Cal¬ 
culations  were  done  for  the  propagation  of  a  shock  wave  through  a  layer 
with  a  semi-infinite  or  finite  thickness  of  dusty  air  containing  10  (im 
glass  particles.  When  the  shock  wave  penetrates  into  the  dusty  gas,  its 
frozen  front  decreases  as  a  result  of  the  absorption  of  energy  by  the  par¬ 
ticles.  Consequently,  the  pressiire  of  the  gas  at  some  distance  behind  the 
shock  front  increases  beyond  its  initial  value  owing  to  the  deceleration 
of  the  flow  caused  by  the  particles.  Similar  transition  effects  occur  for  all 
physical  quantities.  A  shock  wave  is  reflected  from  the  first  contact  region 
into  the  pure  air  and  a  rarefaction  wave  from  the  second  contact  region 
into  the  dusty  gas.  When  the  shock  wave  comes  out  of  the  dusty-gas  layer, 
the  gas  is  released  from  the  restraint  of  having  particles  to  accelerate 
and  the  shock  wave  recovers  much  of  its  strength.  An  expansion  of  the  gas 
occurs  downstream  as  a  result  of  its  acceleration  suid  a  rarefaction  wave 
forms  and  is  reflected  beick  into  the  dusty  air.  The  contact  regions  have 
finite  transition  thicknesses.  For  a  dusty-air  layer  of  finite  thickness, 
the  reflected  rsirefaction  wave  interacts  with  the  first  contact  region  and 
reflects  as  a  conpression  wave,  which  is  transmitted  at  the  second  contact 
region  and  then  overtakes  the  emergent  shock  wave  and  causes  it  to  recover 
more  of  its  original  strength.  For  the  particular  case  of  an  incident  shock 
strength  of  pressure  ratio  of  5.000  the  final  transmitted  wave  strength  Is 
for  all  intents  the  same  (5.005),  Only  the  variations  of  mass  concentrations 
and  temperatures  of  the  mixture  remain  in  the  dusty-air  layer  as  remnants  of 
the  shock-wave  passage.  For  a  weak  incident  shock  wave,  the  transmitted 


shock  wave  in  the  dusty  gas  can  be  fully  dispersed.  Here  the  formation 
of  the  reflected  shock  and  rarefaction  waves  at  the  contact  transmission 
fronts  eire  found  to  take  a  long  time. 

Experimental  verification  of  this  analysis  will  be  done  in  the  near 
firture  at  ITTIAS.  Of  covirse,  different  layers  euid  gases  can  be  used  and 
analysed.  However,  the  present  study  will  provide  inportant  insight  into 
the  forthcoming  experimental  studies  and  into  the  practical  field  cases 
of  blast  waves  passing  throxxgh  dusty-air  layers. 
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PLATE  1(a)  REFRACTION  OF  A  PLANE  SHOCK  WAVE  IN  AIR  THROUGH  A  LAYER 
OF  HELIUM. 


Composite  photograph  of  two  (x,t) -plane  schlieren  records 
showing  the  double  refraction  of  a  plane  shock  wave  at  a 
helium  layer  and  the  resulting  subsequent  overtaking  of  two 
shock  waves  moving  in  the  same  direction. 


/j  SEC. 


PLATE  Kb)  REFRACTION  OF  A  PLANE  SHOCK  WAVE  IN  AIR  THROUGH  A  LAYER 
OF  HELIUM. 

Explanatory  line  diagram  of  (a) .  The  gases  are  at  rest 
and  at  atmospheric  pressure  initially.  The  incident  shock 
wave  Si,  in  air,  refracts  into  helium,  at  the  microfilm 
contact  surface  C2,  causing  a  transmitted  shock  wave  S^, 
and  a  reflected  rarefaction  wave  Ri.  The  transmitted  shock 
wave  St  refracts  into  air  at  the  microfilm  contact  surface 
Cl/  generating  an  emerging  shock  wave  Se,  and  a  reflected 
shock  wave  S,  which  refracts  at  C2,  producing  a  reflected 
shock  wave  Sq;  which  overtakes  Se.  This  gives  rise  to  the 
final  shock  wave  Se,  a  perfect  contact  surface  C3,  and  a 
very  weak  reflected  rarefaction  wave  R2.  The  initial  re¬ 
fraction  at  C2  causes  shock  Si  to  attenuate  in  pressure 
ratio  in  the  helium  layer.  The  second  refraction  at  Ci 
amplifies  the  pressure  ratio  of  Se  in  air.  However,  owing 
to  the  nonlinear  nature  of  the  refractions  Se  is  weaker  than 
Si.  The  overtaking  of  Sq  and  Se  now  amplify  Se  so  that  it 
is  nearly  the  same  strength  as  Si.  Note  the  rapid  diffusion 
of  the  helium-air  contact  surfaces  at  Ci  and  C2  and  the  very 
stable  contact  surface  C3  in  air,  Mach  numbers  of  Si  =  1.70, 
Se  =  1.65  and  Se  =  1.69,  pressure  ratio  across  R2  =  0.97. 


(A)  REFLECTED  SHOCK  WAVE  (B)  REFLECTED  RAREFACTION  WAVE 
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VARIATION  OF  TRANSMITTED  SHOCK  MACH  NUMBER 
Ms  *  Ug/as  WITH  INCIDENT  SHOCK  PRESSURE  RATIO 
P4/P5- 
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FIG.  16  FLOW  QUANTITIES  AT  t  =  0.78xl0'^  SEC  (P4/P5  =  5). 

. —  GAS,  -  PARTICLES. 

(i)  THE  BEGINNING  OF  FORMATION  OF  A  REFLECTED  SHOCK  WAVE 


(d)  Velocity 
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FIG.  18  FLOW  QUANTITIES  AT  t  =  3.12xl0”^  SEC  (P4/P5  »  5) 
■  "  GAS.  -  PARTICLES.  ' 


STATIONARY  STRUCTURE  OF  A  SHOCK  WAVE  IN  A  DUSTY  GAS  (SHOCK 
PRESSURE  RATIO  =  6.09). 
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FIG.  23  -CONTINUED.  , 

FLOW  QUANTITIES  AT  t  »  5.30x10  SEC  (p4/P5  »  1.5). 
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FIG.  24  STATIONARY  STRUCTURE  OF 
(SHOCK  PRESSURE  RATIO  * 
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FIG.  25  FLOW  QUANTITIES  AT  t  =  0.78x10"'*  SEC  (P4/P5  ®  6 

■'  —  GAS,  -  PARTICLES. 

(i)  THE  BEGINNING  OF  FORMATION  OF  A  REFLECTED 
RAREFACTION  V7AVE. 
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FIG.  26  FLOV7  QUANTITIES  AT  t  =  1.56xlo“^  SEC  (p^/p^ 
- -  GAS,  -  PARTICLES. 
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PIG.  30  PLOW  QUANTITIES  AT  t  ■  3.74xl0“^  SEC  \r4/P5  “  1.584). 
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FIG.  30  -  CONTINUED.  - 

FLOW  QUANTITIES  AT  t  =  3.74xl0"'’  SEC  (p^/pj  «  1.584). 

—  GAS,  -  PARTICLES, - EQUILIBRIUM  FLOW. 


CM 


X 

(c)  Temperature 


(d)  Velocity 


.-3 


PIG.  31  -  CONTINUED. 

PLOW  QUANTITIES  AT  t  *  6.86x10  ^  SEC  (p^/pj  =  1.584) 
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FIG.  32  -  CONTINUED.  .  , 

FLOW  QUANTITIES  AT  t  =  2.18x10  ^  SEC  (P4/P5  =  5, 
INITIAL  LAYER-THICKNESS  »  27.2  Cm) . 

-  GAS,  -  PARTICLES. 
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FIG.  33  FLOW  QUANTITIES  AT  t  »  4.99x10  SEC  (p^/p 
INITIAL  LAYER-THICKNESS  *  136  cm) . 
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33  -  CONTINUED 

FLOW  QUANTITIES  AT  t  =  4.99xl0“^  SEC  (P4/P5  “  5, 
INITIAL  LAYER-THICKNESS  =  136  cm) . 

-  GAS,  -  PARTICLES. 


APPENDIX  A 


RANKINE-HUGONiar  AND  ISEMTROPIC  RELATIONS 


We  can  find  the  pressure  P3  =  P2  in  the  region  between  the  transmitted 
shock  wave  and  the  reflected  wave  for  a  given  incident  shock  pressure  ratio 
Plt/P5  solving  the  algebraic  equations  (7)  or  (9).  The  other  flow  quanti¬ 
ties  behind  the  transmitted  shock  wave  can  be  obtained  then  from  the  Rankine- 
Hugoniot  relations  as  follows: 


P2  _  1  +  (7i+l)/(7i-l)  X  (Pg/P^) 
^  (7i^+l-)/(  73^-1)  +  P2^^1 


(Al) 


( 7^+1 )/(  73^-1)  +  Pg/Pj^ 


-  - 

Pl^  1  +  (73^+1)/(73^-1)  X  (Pg/Pj^) 


(A2) 


\  ’'iVPi  /Vv’-/ vPi  ’'i'"^/ 


(A3) 


When  a  shock  wave  is  reflected  at  the  contact  surface,  the  flow  quantities 
behind  the  reflected  shock  wave  are  calculated  in  a  similar  way.  If  the 
reflected  wave  is  a  rarefaction  wave,  the  isentropic  relations  across  the 
rarefaction  wave  are  given  by 


Usii]g  these  relations,  we  can  express  the  values  of  the  flow  quantities  in 
the  \inlform  flow  regions  behind  the  transmitted  shock  wave,  the  contact 
surfaces  and  the  reflected  wave  in  terms  of  the  Incident  shock  pressure 
ratio  or  shock  Mach  number. 


APPENDIX  B 


FORMATION  OF  FULLY-DISPERSED  SHOCK  WAVES 


The  variation  of  the  thickness  of  the  shock  transition  region  with 
time  is  shown  in  Fig.  B1  for  the  case  of  a  fully-dispersed  shock  wave 
(Figs.  21-23).  The  shock  thickness  was  defined  as  the  length  of  the  region 
over  which  the  pressure  varies  from  10%  to  90^  of  the  total  pressure  jun^) 
of  the  stationary  shock  wave  (Fig.  2h).  The  final  thickness  is  also  shown 
by  the  dashed  line.  It  can  be  seen  that  the  trajisition  thickness  increases 
very  gradually.  It  is  estimated  from  Fig.  B1  that  a  stationary  shock  wave 
will  require  a  formation  time  in  excess  of  llxl0"3  sec. 


1 1 


1.56  3.12  4.68  6.24  ^ 

xIO"*  sec 


.  B1  VARIATION  OF  THICKNESS  OF  FULLY 
DISPERSED  SHOCK  WAVE  WITH  TIME. 
( -  FINAL  VALUE) . 
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